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Abstract This paper is concerned with the free boundary value problem for multi¬ 
dimensional Navier-Stokes equations with density-dependent viscosity where the flow 
density vanishes continuously across the free boundary. A local (in time) existence of 
weak solution is established, in particular, the density is positive and the solution is 
regular away from the free boundary. 
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1 Introduction 

The compressible Navier-Stokes equations (CNS) with density-dependent viscosity coeffi¬ 
cients are taken into granted recently. The prototype is the model of viscous Saint-Venat 
system used in geophysical flow [13] to simulate the motion of the surface in shallow water, 
of which the mathematical derivation is also made recently based on the motion of three 
dimensional incompressible viscous fluids on shallow region with free surface condition on 
the top and Navier type boundary condition at bottom of hnite depth [6, 10]. 

In the present paper, we consider the general isentropic compressible Navier-Stokes equa¬ 
tions with density-dependent viscosity coefficients in N = 2,3, can be written for f > 0 

as 

\pt + div(pU) = 0, 

\ {pV)t + div(pU O U) - div(p(p)ro(U)) - V(A(p)divU) + VP(p) = 0, 

where p(x, t), U(x, t) and P{p) = p'^i'j > 1) stand for the fluid density, velocity and pressure, 
respectively, D(U) = |(VU -k VU^) is the stress tensor, and p{p) and A(p) are the Lame 
viscosity coefficients satisfying p{p) > 0 and /i(p) -\- NX{p) > 0 for p > 0. Note here that the 
case 7 = 2 and 61 = 1 in fll.ll) corresponds to the viscous Saint-Venat system. 
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One of mathematical difficulties to investigate the existence and dynamics of solutions to 
fll.ip is that the viscosity coefficients are density-dependent which leads to strong degeneracy 
in the appearance of vacuum [5]. Thus, it is natural and interesting to investigate the influ¬ 
ence of vacuum state on the existence and dynamics of global solutions to fll.ip . One of the 
prototype problems is the time-evolution of the compressible viscous flow of hnite mass ex¬ 
panding into inhnite vacuum. This corresponds to free boundary value problem (FBVP) for 
the compressible Navier-Stokes equations fll.ip for general initial data and variant boundary 
conditions imposed on the free surface. The study is fundamental issue of fluid mechanics 
and has attracted lots of research interests [11, 17]. These free boundary problems have been 
studied with rather abundant results concerned with the existence and dynamics of global 
solution for CNS fll.ip in ID, refer to [4, 9, 14, 15, 19, 20] and references therein. As for 
related phenomena of vacuum vanishing and dynamics of free boundary, the reader can refer 
to [8, 9]. 

The free boundary value problem for fll.ip with stress free boundary condition has been 
investigated in [7], where global existence of spherically symmetric weak solution is shown, 
in particular, the dynamics behaviors and the Lagrangian properties are also established 
therein. Chen-Zhang [3] proved the local solutions of fll.ip with spherically symmetric initial 
data between a solid core and a free boundary connected to a surrounding vacuum state. 
Under certain assumptions that are imposed on the spherically symmetric initial data, which 
between a solid core and a free boundary, Chen-Fang-Zhang established the global existence, 
uniqueness and continuous dependence on initial data of a weak solution in [2]. Wei-Zhang- 
Fang obtained the global existence and uniqueness of the spherically symmetric weak solution 
in [18] with the symmetric center excluded. 

In the present paper, we consider the free boundary value problem for multi-dimensional 
CNS fll.ip in the case that where the fluid density connects with vacuum continuously. 
We show that a spherically symmetric weak solution, with the symmetric center included, 
exists locally in time, in particular the density is positive away from the free boundary but 
vanishes across the initial interface separating fluids and vacuum, and the free surface moves 
as particle pathes in radial direction. To this end, we need to employ the basic energy and the 
modihed Bresch-Desjardins (BD) [1] entropy to establish the expected boundary regularities 
of spherically symmetric solutions in Lagrangian coordinates so as to control the hnite speed 
motion of free boundary within hnite time. Then, in terms of the original equations instead 
of the spherically symmetric form, we are able to apply the higher order energy estimates 
to establish the necessary interior regularities of solutions away from the free boundary but 
with the symmetry center included. Then, the combination of both boundary estimates and 
interior estimates and the above leads to the desired local existence and uniqueness results 
of solutions. 

The rest of this paper is as follows. In Section 2, we state the main results of this paper. 
In Sections 3-5 we establish boundary regularity and interior regularity, with which we can 
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prove the existence and uniqueness in Section 6. 


2 Main results 

For simplicity, the viscosity terms are assumed to satisfy /i(p) = p®, A(p) = pp'(p) — p(p) = 
{9 — l)p^ and D(U) = |(VU + VU^) in fll.ip . The pressure is assumed to be P(p) = p^. In 
this situation, (11.Ih become 

\pt + div(pU) = 0, ^ 

\ (pU)i + div(pU (g) U) - div(p®©(U)) -{9- l)V(p®divU) + Vp^ = 0. 


Consider a spherically symmetric solution (p, U) to (12.11) in so that 

p(x,t) = p(r,f), pU(x,f) = pM(r,t)-, r=|x|, x G (2.2) 

r 

and (12.ip are changed to 

( + {r^~^pu)r = 0 , 

I {r^-^pu), + {r^-^pu\ + r^-\p^)r (2.3) 

[ - r^-\9p%ur + ^u))r + {N- l)r^-\p%u = 0, 

for (r, t) G Qt with 

fir = {(^"P)| 0 < r < a(f), 0 < f < T}. (2.4) 

The initial data is taken as 

(po,Uo)(x) = (po(r),Wo(?^)7), r G [0,ao]. (2.5) 

At the center of symmetry we impose the Dirichlet boundary condition 

M(0,f)=0, (2.6) 

and the free surface dQt moves in radial direction along the “particle path” r = a{t) with 

the stress-free boundary condition 


p{a{t),t) =0, f > 0, 
where a{t) is the free boundary dehned by 

a'{t) = u{a{t),t),t > 0, a(0) = oq. 


(2.7) 


( 2 . 8 ) 


First, we dehne a weak solution to the FBVP fl2.ip - fl2.8p as follows. 
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Definition 2.1. (p, U, a) with p > 0 a. e. is said to be a weak solution to the free surface 
problem fl2.ip - fl2.8p on fit x [0,^], provided that it holds 

0 < p e L°°{t),T-L\flt)nL\flt)), VpU G L^{t),T-L\flt)), 
pWu G L\t),T-L\flt)), a{t) G C\[0,T]), 


and the equations are satisfied in the sense of distributions. Namely, it holds for any t 2 > 
> 0 and (f> G C^{fLt X [0,T]) that 


f p(f)dyi\lttt= [ [ (p0t + pU ■ 

J J t\ J £7^ 


(2.9) 


and for fj = -0^, • • • G C^{flt x [0, T]) satisfying t) = 0 on dflt and V^(x, T) = 0 

that 


mo • '0(x, 0)(ix + 


p'^ div'tpdxdt — (6 — 1) 


p^divlJdivfjdxdt 




'0 JQt 


'0 JUt 


p^VU : Wipd^dt + 


'0 JQt 


'0 JUt 


[pU • dti^ + \/pU (8) \/pU : V'ijj]dxdt = 0 (2.10) 


where mo = mo-. The free boundary condition 1^2. 7]) is satisfied in the sense of continuity. 


Notations: Throughout this paper, C and c denote generic positive constants, Cf^g > 0 
denotes a generic constant which may depend on the sub-index / and p, and Ct > 0 a 
generic constant dependent of T > 0. 

Before stating the main result, we need assume the initial data fl2.5p satisfies for 0 < ro < 
r 2 < ri < rf < oq that 

' jao^jv-i^o(r)dr = 1, 

p*(ao - rY < po(r) < p*(ao - rY, r G [0, Oq], 

' (Po '(Po“Or)r,Po VorWo) e L2([r2,ao]), Uq E {\r 2 , Oq]) , (2.11) 

PoUq™ G L^([r2,ao]), 

^ (po,Uo) G H^{[0,rf]), (y^Uo,Po'^^) G L^([0,ao]), 

where p* and p* are positive constants. 

Meanwhile we list some assumptions on the constants Yi,d,(5,m) with jd = 

(tIi) Let 7 , 6 satisfy 

^ < P < 7, 7 > 1- (2.12) 

{A 2 ) Let ft satisfying 


(2,13) 
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(Aa) m > 0 is a integer satisfying 


m > niax{ 


l + /30-/3’4-4/3 


}• 


Under the above assumptions, we have the following existence result. 


(2.14) 


Theorem 2.1. Let N = 2,3, j > 1. Assume that fl2.1ip and Ai ~ A^ hold, Then, there 
exist a time T* > 0 and > 0 dependent of initial data, so that the FBVP fl2.1l) - fl2.8l) has 
a unique spherically symmetric weak solution for t G [0,T*] 


(p,pU,a)(x,t) = {p{r,t),pu{r,t)^,a{t)), r = |x 


in the sense of Definition lil. il for any T G (0,r*] satisfying that 

pa{t) 


r^ ^p{r,f)dr= I r^ ^pQ{r)dr. 


^ao 


N-1 


(2.15) 


Co < a{t) < 2ao, t G [0,T*], ||a||H2([o,T*]) < C, (2-16) 

(p,U) G X [0,r,]), ||U||H5i.oc(f,,,[0,T*]) < C, (2.17) 

p_{a{t)-ry <p{r,t)<p+{a{t)-ry, (r, t) G [0, a(f)] x [0, T*], (2.18) 

(/9^ +|^Un(x,t)dx+ [ [ p^lVUpdxdt < C, (2.19) 

Jo Jut 

sup ||(p, U)(f)||^3(Qin) + f (||p(t)||^ 3 (Qin) + II VU(t)||^ 3 (™„))(it < U, (2.20) 

ie[o,T.] Jo ^ ‘ ^ ^ ‘ ^ 

pa{t) pT, pa{t) 

sup / pr^~^(pvJ^+ ul)dr + / / p^r^~^{u^^~‘^u^,, + ul^ + r~‘^ul)drdt < C, (2.21) 

iG[0,r*] Jr2;2(0 J^ Jrx2h) 

where fit = {0 < |x| < a{t)}, Q]"' = {0 < |x| < r^^^it)}, r^^f) 'Is the particle path with 
^xi{0) = ri {i = 1,2) and 1 < k < 2m is a integer, and C > 0 is a constant. 


sup 

tG[0,T*] 


'Ut 


Remark 2.1. Theorem \2. 1\ yields the local existence of spherically symmetric weak solutions 
for two/three dimensional compressible Navier-Stokes equation with fluid density connecting 
with vacuum continuously. In particular, it applies to the viscous Saint- Venant model for 
shallow water (which is fll.ll) with N = 2, p{p) = p, A(p) = 0, and ■y = 2). 


3 Basic energy estimates 

The proof of Theorem 12.11 consists of the construction of approximate solutions, the basic 
a-priori estimates, and compactness arguments. We establish the a-priori estimates for any 
solution {p, u, a) to FBVP fl3.4l) - fl3.6l) in this section. 
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Let us introduce the Lagrangian coordinates transform 

p r /*^(^) 


x{r, t) = py dy = l 


py^ ^dy, T = t, 


which translates the domain [0,r] x [0,a(t)] into [0,T] x [0,1] and satishes 


dx 

dr 


= pr 


N-l 


dx 

'dt 


= —pur 


N-l 


dr 


dr 


dr ' dt ’ 


and 


r^(x, t) = N I -{y, T)dy = a{ty^ - N I -{y, T)dy 


N 


lo P 


P 


dr 

d^ 


= u. 


In terms of fl3.ip - fl3.3l) . the free boundary value problem fl2.1l) - fl2.8p is changed to 

pT + p‘^{r^~^u)^ = 0, 

r^~^Ur + (p^ - dp^~^^(r^~^u)^)^ + = 0, 

for (x,t) G [0,1] X [0,r], with the initial data and boundary conditions given by 

(p,M)(a;,0) = (po,Mo)(a^), a:e[0,l], 

u(0, r) = 0, p(l, r) = 0, r G [0, T], 
where r = r(x, r) is defined by 

-^r{x,T) = u{x,t), X G [0,1], t G [0,r], 

and the fixed boundary x = 1 corresponds to the free boundary a(r) = r(l,T) in Eulerian 
form determined by 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


d 

dr 


a(r) = u(l,r), r G [0,T], a(0) = oq. 


(3.8) 


Note that in Lagrange coordinates the condition (12.lip is equivalent to 
p*(l - xY < Po{x) < p*{l - xY, X G [0,1], 
{pY^r^-^UoYx e L‘^{[x2, 1]), pY'^r^-^Uo G iL^([x2,1]), 
urGLi([x2,l]), 

(po,Wo) e if^([0,x|]), {uq,pY~^^^Y G L2([0,1]), 


(3.9) 


where 0 < X 2 = x^~^po{r)dr < xf = JY r^~^po(r)dr. 

First, making use of similar arguments as [7] with modihcations, we can establish the 
following Lemmas l3.1L Lemmas 13.31 which we omit the details. 

Lemma 3.1. Let j > 1, T > 0, and (p, u, a) with p > 0 be the solution to the FBVP fl3.4p - 
fl3.6p forr G [0,T]. Then, it holds 


T /•! 


(b + -1) 

r*r /•! 



p^ ^^dxds 


u 


0 Jo 


+ [1 - N{1 



p^^{r^ ^Uxfdxds = Eq, rG[0,T], 


(3.10) 


0 Jo 
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where 


Eo =: / ilul + ::^pI ^)dx. 


Lemma 3.2. Under the same assumptions as Lemma IS. A it holds 


T /»1 


+ ( 6 * - 1 + ^) / / ^u)^fdxds 


7 

r /•! 


0 JO 


+ 


N-1 

N 



0 JO 


- f^fdxds = Eo, re[0,T]. 


Lemma 3.3. Under the same assumptions as Lemm,a \S. iL it holds 
_ 1 ^ 

iv(7-i) ^iv(^_i) ^ ^ a(r), (x, r) G [0,1] x [0, T], 


7-1 i 


(x 2 - xi)^-i < r"(x 2 ,r) - r"(xi,r), 0 < Xi < X 2 < 1, re [0,T]. 


/n particular, it holds for x = 1 that 


P iV(7-l) 
-^0 


< a(T) = r(l,T), Te|0,r]. 


(3.11) 

(3.12) 

(3.13) 

(3.14) 


Then, we have 


Lemma 3.4. Let T > 0 and 7 > 1. Let {p,u,a) be the solution to FBVP (13.41) - fl3.8p for 
(x, r) e [0,1] X [0, T]. Assume further that it holds for some Xq G (0,1) 

1/9_(1 - x)^ < p(x,r) < 2p+(l - x)^, (x,r) e [xo, 1] X [0,T], (3.15) 

|pr^"^M 3 ;(x,r)| < 2Mo, (x, r) G [xq, 1] x [0,T], (3.16) 

where (3 G (0,1), p+ = 2p*, p_ = Ip* and the constant Mq > 0 is given by L7.,7\) . Then, 
there is a time Ti e (0, T] fl (0,1] so that (p, u, a) satisfies 


cox^(T-i) < r{x,T) < a(r) < 2ao, (x,r) G [0,1] x [0,Ti], 

yt,(x) < p(x,t) < 2fH,(x), (x,t) e [lo.ll X |0,ri], 


(3.17) 

(3.18) 


((/),r"-7'Vs + 



0 Jo 


((p^^')yr^ dyds < CE,„(x,t) e (xo,l) X lO,Ti], (3.19) 


with Ex =: fo(u^ + ^)'^)iy,0)dy + Pq (y)dy and C^o^x > 0 a constant. 

Proof. First of all, it follows directly from fl3.16p and fl3.15p that 


'Xo 


|M|dx+ / \ux\dx 
J Xo 


||'if||L°°([a;o,i]x[o,r]) <(1 ” a:o) ^ 

<Cxo{El'^ + Mop-_^) =:Mi, 

which yields fl3.17p with the help of fl3.12p and 

J'i.x-, t) < o,{t) = Go + [ n(l, s)ds < Oq + TMi < 2ao, 


(3.20) 

(3.21) 

(3.22) 























J. Liu 


for r E [0, Ti^q] with 


It follows from fl3.4p i that 

p{x,t) = po{x)exp 


Ti^a —■ 


{pr 


N-l 


lln 


+ {N — l)ur ){x, s)ds 


(3.23) 


(3.24) 


which together with (13.161) . fl3.17p and fl3.2ip yields (I3.18p for r E [0, Ti^b] with Ti^b determined 
by 

ri,6 =; min{ Ti,,, {2MiE^x~^ + 2Mo)-Hn2}. (3.25) 

Differentiating fl3.4p i with respect to x E [xq, 1), substitnting the resnlted eqnation into 
(l3.4p o and nsing the fact ^ = u, we have 


{u + r^ ^{p^)x)t + ip'^)xr’^ ^ = 0 . 


(3.26) 


Mnltiplying fl3.26p by (p{x){u + r^~^{p^)x), where (j) E C'°°([0,1)), 0 < ^ < 1, (j){y) = 1 for 
y E [0,x] and (piy) = 0 for y > (1 + p)x with y > 0 small enongh, and integrating the 
resulted equation over [0,1] x [0, r] by parts, we obtain 


[ (f){x){u + r^ \p^)^Ydx+ [ [ 

Jo Jo Jo 

<C f (j){x){u +r^~^{p^)^Y{x,0)dx+ 
Jo 

+ {2aof-^M^T{p*f\ 


0(x)((p^2 )^r^ ^Ydxds 



(pxur^ ^p^dxds + 


4>{x)p^ ^{x,0)dx 


(3.27) 


where we have used (13.171) and (13.181) . Choose 


Ti,, =; mini Ti,a, T^^b, {{2aof-^M,{p*f'^)-^Ex}, T, =: min{Ti,,, Ti,,, Ti,e}, 


then the combination of (13.271) and (13.101) yield (I3.19P for r E [0,Ti]. 


□ 


4 Boundary regularities 

This section is devoted to the boundary regularities of solutions to FBVP (I3.4l) - (l3.8p . To 
this end, we first establish the regularities of solution (p, n, a) away from the symmetry center 
and the free boundary. 

Lemma 4.1. Under the assumptions of Lemma \3.4\ there is a time T 2 E (0,Ti] so that the 
solution (p, u, a) satisfies for xi E (xq, 1) and r E [0, T 2 ] that 

IKPx)'U' a;)('?‘)||2,2Q3,|j + f \\{Uxx^Us)(yS)\\i^2(]^^^ rxp^-)ds<C\5^^ (4.1) 

Jo 

\\{Pxxi'^xx^'^t){,'I')\\L'^[[xq,xi\) I II "l^xs) ('S) ||2,2([xo,a:i])^^ — ^5^5^ 

Jo 


(4.2) 
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II {pXXX"} ^XXXl '^XT + / \\{Uxxxx-,Uxxs){ s)\\\'i([xo,xi])ds < CqSI, (4.3) 

Jo 

provided that (po,uo) G H^{I) with [a;o,xi] C (xq where Xq G (0,a;o) and I =: [xq ,x^], 

where Ci > 0, (z = 4,5,6) are constants dependent of Xq and Xi, but independent of Mq, 
64 = ||(po,tio)||Hi( 7 ), ^5 = \\{po,uo)\\h^i), and Sq = \\{po,uo)\\h3{i). In addition, it holds 

\pr^~^Ux{x,T)\ < Mo^a, \u{x,t)\ < Mi^a, (x, t) G [xq, Xi] X [0, Ta], (4.4) 

with Mo,a = { 2 ao)^-^p*{C 4 Sl + C, 6 iy/^ + Mo,aPf^). 


Proof It is easy to verify that fl4.4p follows from fl4.1l) , fl4.2l) , fl3.18p and (13.201) . What left is 
to show dHP-dMl). Rewrite ([MD 2 as 

r^-^Ur + (p" - 9p^^\^-\x)x + il-9)iN- l){p%- - eiN - l)p%-)x = 0, (4.5) 


Take inner product between (14.5p and (j)p^~^Ur, where cj) = 'if‘^{x) and if G C'°°([0,1]) satishes 
0 < "0 < 1, "0 = 1 for X G [(1 — 2?7 )xo, (1 + 2p)xi], and -0 = 0 for x G [0, (1 — 3?7)xo] U [(1 + 
3p)xi, 1] with a hxed constant p G (0,1) small enough so that [(1 — 3p)xo, (1 + 3?7)xi] C 
(X(C,x0). By lemma [331 and a direct computations, it follows 


i Jo ^Ux)dx + ^ (j)p^ ^r^ ^uldx 

< Cx,{Mo + Ml + {py-^) y'^+^-^r^-^ldx + ao,xiMoMi(p, + p*)i+^-^ 

Jo 

P Cxo,xia^ ^(p* + P*)^^"*"^ ^ + C^xo.xi^o (P* + ^ P Cxf,,xiMi{p* P P* 

+ Cx,,xtPl~''ExtiM^O P Ml P (p. + p*)27+l-iV-2.)^ 


2e-i-v 

(4.6) 


where C^o.^i > 0 is a generic constant dependent of Xq, Xi, but independent of Mq, Mi. 
Integrating (14.6p over [0, r] to get 




'0 



0p^ ^u^dxds 


0 JO 


(fp 


<Cx,,xAMoPM, + {pyy 

>0 JO 

A-N T7 t T\/r2 I T\/r2\\ , ^2 



1 


2+0-V.„V-l„,2 


^wfdxds + t{Cxo,xAMo + Mf + Oq ^) 


+ C'xo.xiP* ^xfiMo + M;^)) + Cxo,xi54, 


(4.7) 


which together with Gronwall’s inequality leads to 


f ^r^ ^uldx + f [ 0p^ ^r^ ^u^dxds < Cxo,xiA^I- 

Jo Jo Jo 

Due to the fact that 0(x) = 1 for x G [(1 — 2p)xo, (1 + 2p)xi], it follows 

i>{l+2ri)xi i-T i‘{l+2ri)x\ 

/ uidx + / / utdxds < C'xn^, 25l, 


(4.8) 


(4.9) 
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for some constant C 3 ;g^xi ,2 > 0 and r G [0, T 2 ^a\ with T 2 ^a chosen as 

T 2 ,a =: min{ri, ln2 }, 

Ah =: + Ml + a^-i) + + M^), 

K 2 =-.C,,^,,{Mo + Ivh + {py-<>). 


In addition, it follows from (14.Oh . (Id.dh o and Lemma [3.41 that 

n (l+2»7)xi 
,L-2ry)xo 

The combination of fl4.9l) - fl4.10l) and fl3.19p gives rise to 


(^xx A — C'xo,xi<54. 


(4.10) 


^(1+2?7)xi i>{l+2r])xi 

/ {ul + pl)dx+ / {ul^ +pl^ + ul)dxds 

' {l—2ri)xo Jo J (l—2r))xo 


(4.11) 


for r e [0,T2^a], which implies fl4.ip for any T 2 < T 2 ^a- 

The higher order regularities of the solution can be obtained by applying the similar 
arguments as the proof of fl4.11l) . Indeed, differentiating fl4.5p with respect to r gives 


^Urr — (1 — N)r ^'^UUr + {p^ — Ux) 


^-N 


e+iN^ 


u 


+ {l-e){N- !)((/).-). - 9 {N - l){p%^)x)r = 0 


(4.12) 


taking inner product between (I4.12p and (pUr over [0,1], where (p = xp‘^{x) and rp G C'°°([0,1]) 
satishes 0 < 'ip{x) < 1, ip^x) = 1 for x G [(1 — p)xo, (1 + p)xi\, and xp{x) = 0 for x G 
[0, (1 — 2p)xq] U [(1 + 2p)xi, 1], we can obtain 

<pr^~^uldx + ^ (pp^^^r^-^ul^dx 

= [ (px{p^ - 6p^^^r^~^Ux)rUrdx ^ f (pr~^uuldx + 6{N - 1) f <p{p^ {^)x)rUrdx 

Jo 2 Jo Jo 

+ [ (P{{p^)t - 9{p^"^^r^~^)rUx)uxrdx -^ [ (pp^^^r^~^ul^dx 


>0 


>0 


+ {l- 9 )il-N) / (P{{p%^)xUxdx 


^Cxo,xi{^ Ml-\-(p^p*Y ^) / (pr^ ^updx -\- CxQ^xii^ {p*Y~^^ 

Jo 

r{l+ 2 r])xi x(l+ 2 r!)xi 

+ {p*+P*Y ^) / uldx + Cxo,xi{p* + P*Y ^ / plrdx 

J {l— 2 ri)xo J {l— 2 rf}xo 

+ ao,x4,3(i + (p*)'^"+'^ + (p*)')(i + K + Ml), 

which together with Lemma [3.41 and 04.111) yields 

f (pr^~^uldx+ f [ (pp^^^r^~^ulJxds < Cxo,xidl, rG[0,r2,d], 

Jo Jo Jo 


(4.13) 


(4.14) 
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with T 2 ^b chosen as 

T2,b =■■ mm{T2,a, K^^Sl + M, + {p, + p*f-^))-Hn2}, 

=: a.....,s(l + (p*)2(T'+h + (p*)4)(l + + Mf). 


(4.15) 


^3 • ~r \P j 

This and flS-dp o. fl3.17p - fl3.18p imply that 

|•{l+r])xl rr p{l+ri)xi 

/ iuL + ul)dx+ / ulJxds < re[0,T2,fc]. 

J {l—ri)xo Jo J {l—r])xo 

Meanwhile, taking inner product between fl3.26p .^ and (f){x){u + r^~^{p%x)x over [0,1] x 
[0, r], and using Lemma [131 we can obtain 

^ (j){u + r^-\p^)^)ldx+^ (t)p^-\u + r^-^{p^):,)ldxds 

= j (f){u + r^~\p^)^)l{x,0)dx+‘^ [ [ (l)p'^~\u + r^~^{p^)^)^u^dxds 

r*T /•! 


0 JO 


6 



(pr^ ^{u + r^ ^{p^)x)x{p'^ ^)x{p^)xdxds 


0 JO 


(4.16) 


+ (p. + p*)’-’ + (p. + p*)"'"-'’ 

+ (p*y^'’ [ max[(p^)^r^“^]^ds)^|. 

Jo 

Integrating fl3.26p over [0, r] to get 

’’^"^(P^)x(a;,'r) = r^"^(a;,0)(po)x - [ {p'^)xr^~^{x, s)ds - u{x,t) + u{x,0), (4.17) 


one deduces from (13.181) . fl3.2ip and fl4.16p that for any (x,r) G [(1 — 2p)xo, (1 + ‘2p)xi] x 

[0, T2,b]-, 

[ [{p'^)xr^~^?{x,s)ds 
Jo 

^*^ 2 {-y-e) f f s)dzds, (4-18) 



0 Jo 


which implies for (x, r) G [(1 — 2p)xQ, (1 + 2?7)xi] x [0, T 2 x] that 

[ [{p'^)xr^~^T{x, s)ds <2ln2, rG[ 0 ,T 2 ,c], 

Jo 

with T 2 x chosen as 

T 2 ,, =: min{ r 2 , 6 , {C{p*f^'^-^\l + M^))-Hn2, {C{p*f^'^-^^)-Hn2}. 
Using Lemma 13.11 fl3.19l) and fl4.ip that 

p(l+r;)a;i pt p(1+T7)xi 

/ Plxdx+ / plJxds < Cxo,x,dl rG[ 0 ,r 2 ,c]. 

J {l-r))xo Jo J {l-ri)xo 


(4.19) 


(4.20) 
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The combination of fl4.15p . fl4.20p . fl4.ip and fl3.4l) o gives rise to 


•>(l+77)3;i 


r i'{l+r])xi 


'{i--n)xo 


(«r + “xx + Pxx)dx + 



(“Lx + pIx + ul,)dxds < (4-21) 


0 J {l—ri)xo 


for r e [ 0 ,T 2 ^c], which implies 04.21) for any T 2 < T 2 ^c- 
Differentiating 04.51) with respect to x to get 


UxxxT + ^u^)x) 


+ (AT - 1)(1 - e){r^'^-\pXuU^ - e{N - l)K-VLy)x)xxx = 0, (4.22) 

and taking inner product between 04.22p and (t){x)uxxx over [0,1] x [0,r], where 0 = 
and ?/) G C^([0,1]) satishes 0 < i/j^x) < 1, i/j^x) = 1 for x G [xo,xi], and 'i/j{x) = 0 for 
X G [0, (1 — p)xo] U [(1 + p)xi, 1]. We can obtain 


1 d 

2 dr 


(pulxxdx- / (f){r^ ^U^)x)xxUxxxxdx 


= - ^(p'^-p^+V^ ^U^)x)xUxxxxdx - / (pxx{r^ ^(p'^-p^+V^ ^U^)x)xUxxxdx 


+ 9{N-1) I (t){r^ ^p%^)x)xxxUxxxdx + {N - 1){1 - 6) / (j)xUxxx{r^ '\p^)xu)xxdx 


+ (A^ - 1)(1 - 0) / (j)Uxxxx{r \p^)xU)xxdx. 

Jo 

Lemma 13.41 03.4p , 04.111) and 04.2ip lead to that 

^ I (t){ulxx + lP^~^plxx)dx+ I (t)p^^^r‘^^^-^\lxxxdx 


(4.23) 


C*X 0 ,X 1 (1 + <^4 + <^5 + ||U2;3;(t) p_|_^)3,j])) / 4’{p'^ Pxxx '^xxx)dx 

Jo 

+ Cxo ,xi(l T ^4 T ^ 5 ) T C*xo,xi ||'lfxx('^) ||l°°([(1— r;)xo,(l+» 7 )xi])h 5 , 
we apply the Gronwall’s inequality, and the fact 'i/j{x) = 1 for x G [xo,xi], we can obtain 


r-xi 


(“xxx J~ Pxxx')dx T 


T /-xi 


'xo 



ulxxxdxds < Cx^,xA^ rG[ 0 ,T 2 ,d], (4.24) 


0 J XQ 


with T 24 determined by 

^ 2 ,d =: minj T’ 2 ^c) (C*xo,xi(1 + ^4 + (C*xo,xi(l + <^4 + hg)) ^ln 2 }. 


By fl434l) and (I 33 P 2 , we have 


PXi 


T nx\ 


ul^dx + 


fXQ 



ulxJxds < Cxo,xi5l, T G [0, r 2 ,d]. 


0 JXQ 


The combination of 04.241) and 04.251) yields 04.3p with T 2 = T 2 ^d- 


(4.25) 

□ 
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Lemma 4.2. Under the same assumptions as Lemma 3.4. there is a time T 3 G ( 0 ,T 2 ] so 
that it holds for (x 2 ,t) G (xo,Xi) x [OjTs] 


T /»! 


T)dx + 


' 2:2 



pe+ 1 ^ 2 k- 2 ^ 2 N- 2 ^ 2 j^^^ < Cy, = 1 , 2 , • • • , 2 m, 


0 J X2 


T n\ 


+ p® ^u^)dx + Ugdxds < Cs+ p^^'^u'^^dxds, 



T /•! 


0 ./X 2 





<C, + C{I (p^+^ui + p^-^u^ + ui)dxy, 


0 J X2 


fi+1^,2 I ^6J-K,2 I ^ 2> 


^X 2 


T pi 


'X 2 


e+ 12 N- 2^2 , 2^ 


uidx+ / / {p^+W^^-\i, + p'^-^r-X)dxds<Cio, 


(4.26) 

(4.27) 

(4.28) 

(4.29) 


0 J X2 


where Ct = Qdlpollnid^o,!]), ll^io||i^ 2 ([xo,i])) > 0 are constants, i = 7,8,9,10. 

Proof. We apply the arguments used in [3] to show 04.261) . First we consider the case of 
k = 1. From Lemma 13.11 we obtain 04.261) with k = 1 easily. Assume 04.26P holds for 
k = I — 1, 


T pi 


' X2 



“^uldxds < C. 


(4.30) 


0 J X2 


Now we need to prove 04.261) holds for k = 1. Multiplying 03.41) 9 by and integrating 

over X from 0 to 1, where (j) = '^^(x) and '0 G C'°°([0,1]) satishes O<'0^1)'0=1 for 
X G [x 2 ,1], and V’ = 0 for x G [0, Xq]. We can obtain 

d 1 


(fu^^dx 


dr Jq 21 
-1 


= — 0 (j)p^~^^{r^ ^)x{r^ ^u)xdx+ / (f)p^{r^ ^u^’' ^)xdx 


'0 


+ {N-1) cPp%P^-%^%dx-e p'^+V^"-V‘-Vx(fo'“'«)xdx 

Jo Jo 

rl /•! 


+ J p'^r^ ^(pxdx + {N — 1) I p^r^ ‘^u^^'c^^dx 
:=Hi + H2 + H3 + H^ + H5 + H^. 


'0 


(4.31) 


Set 


thus 


Bl = p"+'u2Z-2^2V-2^2 > > 0, 


6'-l„-2„,2i 


+H 3 = -e{2l-l) (j)Bldx-2{N-l){e-l)l (j)BiB 2 dx 


+ {{N-l){N-2)-e{N-lf) / c^Bldx. 

Jo 

Inserting this to 04.311) and using Young’s inequality, we get 


A ri 

dr In 21 


(fu^mx + 9if2l — 1) / cfBfdx 
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<e / (l)Bldx + Ce / (j)Bldx+ / (j)p^{r^ ^):^dx + + \He\. 


By a direct computation, it follows 


J (pu^^dx + J ‘^ul.dx 


<C^,{1 + + Mf-^) + a„(l + Ml) I cPu^^dx 


2l-2\ 


,2i. 


using Gronwall’s inequality and the fact that 0 = 1 for x G [ 3 ^ 2 ,1], we get 

f u‘^^dx+ [ [ p^~^^v?‘^~‘^r‘^^~‘^ul.dxds < C, rG[0,T3^a], 

J X2 J 0 J X2 

with T 3 a determined by 

T 3 ,a =; min{(a„(l + + Ml^-^Mo,a + (^^(1 + Ml))-Hn2} 



we get (14.2611 immediately. 

To show fl4.27p we multiplying (13.411 9 by (pr^~^Ur, where (p = 'ip'^{x) and ^p G (7°^ 
satisfies 0 < -0 < 1, "i/; = 1 for a: G [x 2 ,l], and '0 = 0 for x G [0,Xo], integrate 
[0,1] X [0, r] to have 

T pi 

<pul{x, s)dxds 

0 30 

r /»! nr n\ 

(pp"*{usr^~^)xdxds — 9 / / (pp^^^{r^~^u)x{usr^~^)xdxds 

0 30 Jo Jo 

n l PT p1 

(pp^{r^~‘^uus)xdxds + / (pxp"*Usr^~^dxds 

Jo Jo 

n l PT P 1 

(pxp^^"^{r^~^u)xUsr^~^dxds + {N — 1) / / (pxP^r^~‘^uUsdxds 

Jo Jo 

T pi pi 

6>+l, 2 I 0-1^2\ 



<CSl + C I / (pp^^^ul.dxds — Cl (p{p^~^'^ul. + p^ ^u^)dx^ 
'0 Jo Jo 


which implies 



T pi p1 

ul{x, s)dxds + / + p^~^u^)dx 

0 J X2 J X2 

rr i‘X 2 rl 

<C + C { + )(Pp^^\tdxds 

Jo 3x0 J X2 

<C + C f [ p^~^^u‘l.dxds, 

Jo 3x2 


last we get (14.2711 , where we use the Lemma 13.21 Lemma 14.11 (13.1911 and m > 


(4.32) 


(4.33) 


>([ 0 , 1 ]) 
it over 


(4.34) 


9)- 















FBVP for CNS 


15 


Now we begin to prove fl4.28l) . Multiplying fl3.4P o with 0, where (j) = and '0 G 

C'°°([0,1]) satisfies 0 < 0 < 1, 0 = 1 for x G [x 2 ,l], and 0 = 0 for x G [0,Xo], then 
integrating over [x, 1 ], we have 

(j)Ux — {N — l)(j)p~^ur~^ + — l)0p“^nr“'^ 

u u 


1 

1 


+ l(iV-i)p-i-V-^ / 0(^)ydy - / cPr^-^Uxdy 


1 

e' 


-1 


+ j - ep^^\r^-\)y)dy 

+ f pUy^dy, 


integrating it over [ 0 , 1 ], we have 


^ 61 + 304 ^ 4 ^^ 


< C / 0 '^p^'’' ^dx + C / (j)‘^p^ ^u'^dx + C / p^^^i (j){—)yP^dy)^dx 


,u. 


>0 Jo 

ri rl 

,-361-1/ / „1-N 


pel p ^''0ns(i|/)'^dx + (7 / p^^ ^{j p^tpy—dy^dx 

Jo J X Jo J X ^ 

+ C [ p~^^~^{ f (j)y{p'^ - 6 p^^^{r^~^u)y)dy)'^dx 


u 


'0 


1 p1 p1 

2X+i^,2 T \2 I r^r I j,2 „e-l,,2^„^2 , / j,2,, ^,,^2 


<C + C{ ^^p^+WjxY + C{ / (P^p^-^u^dxY + C'( / 0^n,di/) 

Jo Jo Jo 


which implies 


and 


[ (p'^p^'^^u^dx < C + C{ f (Y{p^ + ul + p^~^^ul)dxY, 

'0 Jo 


f P<>^\ld <C + C{ [\p^-^u^ + ul + p^^WYdxY. 

JX2 JX2 

Finally we prove 04.291) . Differentiating equation 03.4l) o respect to r, 

Urr = {r^~^{6p^"^^{r^~^u)x - p'^)x)t -{N - l){r^~‘^u{p^)x)r = I + J, 


(4.35) 


(4.36) 


/ = r^-\ep<^^\r^-\Y)^ -{N- l)r^-‘^uYpY 


_ (AT _ l)(iV - 2)r^-\\p^)x -{N- l)r^-Mp').r. 


where 
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Multiplying fl4.36p by (pUr and integrating over x from 0 to 1, where (p = pP{x) and 
^p G (^“([0,1]) satisfies 0 < ?/> < 1, ?/> = 1 for x G [ 2 : 2 ,1], and '^ = 0 for x G [0, Xq], we have 

d (p 


1 

s+i\r^N-i„, \ i2 I N-i / 


-uldx + + 6 - 1) / ^Ur)x]‘^dx + 


1 I r) “T - - \ N ' I-''' ^Tjx\ I JV I P dx 

dr Jo 2 Jo Jo 


= — f 9(pxP^~^^r^ ^Ur{r^ ^Ur)xdx + {N — 1) f p^CpxT^ ‘^U^dx + f CpJUrdx 

Jo Jo Jo 

<£ / (pp^~^^r‘^^^~^'*ul^dx + e ( (pp^~^ulr~‘^dx + Cs f (pp^~^^u‘^dx 

Jo Jo Jo 

+ C [ (Pu'^^dx + C f + C /" (Pp^^^-^^+C, 

Jo Jo Jo 

after integrating over [ 0 ,r] and choosing proper e, which implies 

f uj-dx + [ [ + P^~^r~‘^u 1 )dxds 

J X2 Jo Jx2 

n l pi 

p^~^^u‘^dxds + / ul{x,0)dx + C 

2 J XO 

<C [ { f {p^~^vJ + Mg + p^^^ulJjdxYds + C, 

Jo J X2 


Using (14.271) and above we have 


(\p^^w^+p^-w+ui)dx+ r +/-v-x) 

X 2 Jo J X 2 

T pi 


< Uii + Cu (/ + ui)dxyds, 

Jo J X2 

by Gronwall’s inequality, we get 

pi pT pi 

^0+1^2 I J-1,.2 I ,,2\i^ I / / / 6>+1^27V-2, 2 , ^0-1 -2, 2 


^a :2 

where 


(p'^+X + P^V + <)dx+/ / + 0<r<T3, 

0 Jx 2 


Tg =: min{r2,rj 


1 


3,0) 


2 C'iiC'i 2 


-}, 


thus we get (I4.29p . 


□ 


Lemma 4.3. Under the same assumptions as Lemma \3.4\ the solution {p,u,a) satisfies for 
the X 2 G (xo,Xi) and r G [OjTs] that 

[ \uxy°dx < C, (4.37) 


'X 2 


where Aq is a constant satisfying: 


1 < Ao < min{ 


4m 


4m/5 + 1 ’ fi{6 + 1 ) 


}• 


(4.38) 
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In addition, it holds 

I |< Mo, 6 , \ u{x,t) \< Mi^b, (x, r) G [xa, 1 ] x [ 0 , T 3 ], 

where Mq^i, and Mi^b are given by l\4-46\ ) and U.44\ ) respectively, and 

peC'°([x 2 ,l] X [0,T,]), ueC%[x2,l] X [0,T3]). 

Proof. From (13.41) 0 and boundary conditions, we have 


(4.39) 


(4.40) 




=-r^ ^p"' ^ ^ — {N — l)p ^ur ^ + ^{N — l)p ^ur ^ 


+ -{N -l)p-^-^r^-^ 


e 


r^ ^Urdy, 


(4.41) 


u„ 


from fl4.41D . we obtain 

,\^°dx<C [ {l-x)^°^^'^-^-^^dx + C{ f p~^^^°dx)^^{ f u^^dx)^ 

J X2 J X2 

(1 - 

+ C(/ Z+^dp)^ f \p-^-\f p^-Hy)-d\^^dx 

J X2 •J X2 'J X 

+ C{[ p^-\^dy)^ [ \p-^-\ [ p^-^dy)^\^°dx <C: 


+ C'||Wr||i2[j,2^i] 


'13, 


(4.42) 


'X 2 


’X 2 


SO we obtain fl4.37p . where we use the fact Aq < min{ ^^g_^^^ , 4 ^^^ }. Then we have 

1 


\U\\L''°{[x2,1]x[0,T3]) < Y 


X 2 


|ri|dx 


\ux\dx 


^ Cx2 Eq 


1^x1^° dx) ^0 [ I dx) ^0 
J X 2 

< Cx,eI + Cx,C^° =: Ml,6, 


' X 2 


with 


Mi,b =: maxiCx^E^ + Cx^Cf^ , Mi,^}. 
From fl4.4ip for x G [xa, 1] we have 

\pr^~^Ux\ =\\p"'~^ — {N — l)ur~^ — l)ur~^ 

o u 


(4.43) 

(4.44) 


+ t(Ar-i)p ® I ^l)yp0dy-^p ^ I ^Urdy\ 


,u, 

r' 


1 


<C + C{py-^ + C{1 - x)5-®^ + C{1 - a;)|(/3('?-i)+i)-'?/3 

< Mq 6, 


(4.45) 
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with 

Mo,b =■■ max{C + C{p*y-^, Mq, J, 

which and (14.431) lead to (I4.39p . 

We turn to prove fl4.40p . It is easy to verify 

e L-(0,Ts, H\[x2, 1])), ip^)r e L“(0,Ts, L^{[x2, 1])). 


(4.46) 


(4.47) 


This implies p'^ E C'°([x 2 ,l] x [OjTa]) and the continuity of density p E C'°([a: 2 ,l] x [OjTs]). 
Indeed, it follows from fl3.4p i 

(p^)r = = -7P^+V^-i«,-7(iV-l)p^wr-' G L-(0, Tg, L2 ([x2, 1])). (4.48) 

On the other hand, one derives from (14.17p 


'X 2 


{r^-\p^).fdx<C / {p^-^p,^fdx + C{p*f^^-^^ 

J X2 

+ C{p*f^'^-^^ f iu^ + ul)dx 

J X2 

< (7(p*)2(7-e) + C{p*f^'^-^^Eo 


T n\ 



(r^ ^{p"*)xYdxds 


0 J X2 


T /»! 


+ C{p*f^'^-^n {r^-\p^)^fdxds, 


(4.49) 


0 J X2 


and then 


{r^~\p'^)x?dx < C(p*)2(t'-®)(Eo + 1)(1 + '^), (4.50) 


•X 2 


which implies 


[{p'^)xfdx < C'(p*)2('>-®)(Eo + 1)(1 + '’), (4.51) 


fX 2 


thus we obtain 

p^eL-(0,T3,i7i([x2,l])), 

this and fl4.48p gives the half of fl4.40l) . We can also show 

U E L°°(0, Tg, W^’^{[X 2 , 1])), Ur E L°°(0, Tg, L2 ([x2, 1])), 

for any p e fhaf 


(4.52) 


(4.53) 


sup \\u4lp{[x2,i]) < sup ||(pr ) ^LH[x 2 ,i])-\\pr (4.54) 

re[o,r3] te[o,T4] ^ ‘ 


and ^ 

\\Ur{'r)\\l2[a:,,l]= uldx<C, 

J X 2 

this and fl4.54p implies the continuity of velocity u on [x 2 ,1] x [0,Tg]. 


(4.55) 

□ 
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5 Interior regularities 

It is convenient to make use of fl2.ip directly to investigate the interior regularities of solu¬ 
tions. Indeed, we have 

Lemma 5.1. Under the assumptions of Theorem \2.1[ there is a time G (OjTa] so that the 
solution (p, f/)(x, f) = (p(r), ti(r, t)|^) to the FBVP H2.1\) and H2.5\) satisfies 



where =: {0 < |x| < rx 2 ii)}, is the particle path with rx2{0) = r 2 G (ro,ri), and 

M > 0 is a constant given by fl5.14l) . In particular, it holds 



(5.2) 

(5.3) 


with Cs> 0 the Sobolev constant for ||/||l°o < C's||/||_h -2 . 


Proof. We hrst choose < T 3 to be small and assume that it holds 



with 


A/„ = CoM'/\ 


(5.5) 


and M 2 = . It follows from fl5.4l) and (13.171) 


M(t)||L°o([0,r^2(b]) “ l|U(t)||l,°°([0,r^2(t)]) < 2ao || VU(f) ||iooWl) “ ( 5 . 6 ) 


which together with (14.4p . (14.391) . (13.181) and (13.241) (or (I2.1l) ii yields 


\u{r,t)\ < Ml, (r,f) G [0,a(f)] x [0,Ti,a], 
^Po{r) < p{r,t) < 2po(r), (r, f) G [0,a(f)] x [0,r4,J 


(5.7) 

(5.8) 


with Ml =; max{Mi_b, 4aoMo}. 


Corresponding to (13.101) . we have the basic energy estimates 



(p |Up-I-p'^)dx-I- / / (|p2 VUp-I-|p2divUp)dx(is < i?( 


' 0 - 


(5.9) 


0 Jo Jo 


Take derivative 5“ with 1 < |q(| < 3 to (12.ip to get 


(9“log/3), + U . Va“logp + V . S“U(x((),() = 
9"U, + U . va“u + 7 p^-‘V9“ log p = + /„ + 


(5.10) 

(5.11) 
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where 

ga = -<9“(U • Vlogp) + U • Va“logp, K = a“(p-Miv(pVU) + {e- l)p-V(p®divU)), 

/„ = -(a“(u-vu)- u-va“u), k = -7(5“(p^”Viogp-p^~Va“iogp). 

Take inner prodnct of 05.111) and (j){x{r,t))d°‘\J over [0,a(f)] x [0,f], where (f) = and 

G (^^([O,!]) satishes 0 < ?/>(?/) < 1, 'ijj{y) = 1 for p G [ 0 ,X 2 ], and i/j{y) = 0 for p G 
[(1 — p)xi,l] with g > 0 small enongh so that (1 — g)xi > X 2 , and nse the facts that 
4>t = (pr = 4>xpr^~^, make nse of 05.101) and the relation 

/ logp =\jg. J logpp — I + V • (0p'’'“^U)]|(9" logpp 

— f (j)gap'^~^d°' logp — f d°‘ logp5"U ■ 

we have 

pap) 1 

/ ^p^-^\d‘^ log p\'^)r^-^dr 

Jo ^ 2 

/•“(o) 1 ^ 

= / cP{-\d^\J\^ + ^p"^-^\d^logp\y^-\r,0)dr 

In ^ ^ 


nt ra(s) 



(j)s\d’^V\‘^r^-^drds 


t /•<2.(s) 


'0 JO 
(•t /•^'(-s) 



‘U ■ (U ■ Va“U)r^-Mrds 


0 Jo 



'0 JO 

ft l‘a(s) 



[{(pp'^ ^)s + V-(^p'^ ^U)]|5“ log ppr'^ ^drds 
d°‘ logp(9"U • V{(j)p'^~^)r^~^drds 


0 JO 
t r(i(s) 


+ 



0 JO 


09“U(/tQ, +/q, + fcQ)r^ ^drds 

t /•a(s) 

(pr^~^gaP^~^d°‘ log pdrds, 



0 JO 


after a direct compntation that 

pap) 

/ (/)(|a“U|2 + p^-i|a“logpnr^-idr 


t /*a(s) 


+ 



0p®“^(|Va“U|2 + |V • d^V\y^''-^drds 


iaTT|2\„Af-l, 


0 JO 


<C||(Uo,Po)||^ 3 ([o,ri]) + Cxo,xi{l- + Mq + Ml + M 2 + {p*y ^)((54 + ^5 + hg)f 

+ Cx,,xAM^ + Ml + M^ + 52 + ^ 2 ) r 0(|a“U|2 + p^-i|a“ logpHr^-idrds 


'0 JO 
pt 3 /*a(s) 


+ ^7(1 + Mo + M 2 + 54 + hs) / (5^ / 0(|a“U|2 + p'>-i|a“logpnr'^-Mr)2ds, 

l«l=i' " 


'0 


'0 


which implies 


^(0 +J2f f 0p®”^(| Va“U|2 + |V ■ <9“unr^-idrds 


|a|=l 


0 Jo 
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<Ke + K7t + Ks / Y\s)ds, 


(5.12) 


where 

Kq = C||(Uo,Po)||H3([0,ri])> 

+ M 2 + ^4 + (Jg) + (7(1 + Mq + M 2 + ^4 + hs), 

K? = CxQ^xii^ + Md + + M 2 + (p*)"'” ^)(h4 + hg + (5g) + Cxq,xi{Mq + Mf + M 2 + ^4 + (Jg), 

^4,b =: min{i7f^(7||(Uo,Po)||^3([o,^j]), T^^a}, 

and 


-.a(t) 


Y{t)= / (/)(|(9U,(9^U,9'^U)f + K^logp, 5"^logp, 5'^logpf)(r, ^dr. 

Jo 

We apply Gronwall’s inequality to have 

/■“(*) 

0(|<9“U|2+la^logpHr^-^dr 

|a|=M0 


JL Y /■“(*) 

_j_ \ ^ / / J. l/'IV70aTT|2 

|a| = l'^0 io 
<4(7||(Uo,Po)||^3([o,ri])) ^ ^ [0,714 




(5.13) 


for T 4 =: minj r 4 f,, r 4 ^ } with T 4 ,, = (di^g-^s) The fl5.13p together with fl5.8l) . Lemma l4lT] 
and 

lldllH3(0™) — ^3;o,Xi (II logpll j^3(Qin)) , 

leads to fl5.ll) with M given by 

M = ao,.i(4C||(Uo,po)||^3(^..))^ (5.14) 

The estimates fl5.4l) and fl5.2l) follow respectively from fl5.1l) and the Sobolev embedding 
theorem, and fl5.3p follows from fl3.24p and fl5.ll) . □ 


6 Proof of the main results 

Proposition 6.1 (Existence and Uniqueness). Under the assumptions of Theorem, \2.1\. 
there exists a time T* > 0 dependent of initial data, so that the FB VP (12.ip and (12.51) admits 
a unique solution 

(p,pU,a)(x,t) = (p(r,t),pn(r,t)^,a(t)), r = |x|, {r,t) e [0,a(t)] x [0,T*], 
which satisfies 

IIVU||ioo([0,r^2(p]x[0,r*]) < Mq, \\pr^ ^'*^a:||L°°([a;o,l]x[0,T,]) < Mq, ( 6 . 1 ) 
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Ipo{x) < p{x,t) < 2po{x), (x,t) e [0,1] X [0,T*], 


( 6 . 2 ) 


for constant Mq > 0 dependent of initial data, p+ = 2p* and p- = In addition, 

(p, pU, a) satisfies the inner regularities in Euler coordinates 


fa{t) 


(p|U|2 + p^)rfx + 

3 


t 



p^\VUfdxds <C, t E [0,T*], 


0 JO 


|a| = l 



(|a“pp + < ^n, t E [o,r,], 


0 Jni 


(6.3) 

(6.4) 


where LVf" = {0 < |x| < r^^it)}, r^^iit) is defined by r(,^(t) = u{rx^{t),t) with ^^.^(O) = ri E 
(r 2 , ao) and xi = 1 — por^~^dr E (x 2 , 1), and C > 0 and Cin > 0 are constants, and the 
boundary regularities in Lagrange coordinates 


T /»! 


J+1^2k-2JlN-2^2 


u^^^dxE / p^+V-^r^^-^uidxds<Cb, re[0,T*], 


0 J X2 


r pi 


0+12N-2^,2 


uidx + / / + p'^-^r-X)dxds < a, re [0, T,], 


0 J X2 


(6.5) 

( 6 . 6 ) 


with the integer 1 < k < 2 m, and Cb > 0 is a constant. 


Proof of the Theorem \2.1\ With the estimates we have obtained in Sections 3-5, we can 
apply the method of difference scheme and compactness arguments as in [3, 8] and references 
therein, to prove the existence of weak solutions to the FBVP fl2.ip and fl2.5p . we omit here. 
Next, we apply the idea in [3] to prove the uniqueness. Let (pi,rii,ri) and {p 2 ,U 2 ,r 2 ) are 
two solutions to the FBVP fl3.4p - (l3.8l) . and denote 


Ti 

{g,uj,R){x,r) = (pi - p 2 ,tii - U 2 , -l)(a:,r). 

r2 

Based on Proposition 16.11 we can derive the following estimates 
1 pi(x,r) P 2 (x,t) 


0 < cp*(p*) < 


+ 


< C'p*(pJ \ (x,r) e [0,1] X [0,T*], 


P 2 {x,t) pi(x,r) 

— — -h \p\'^^r^~^Ui,,\ + |p 2 +Vf"^M 2 x| <c, x E [0,1], 

ri r2 


/ ri(x,r) r 2 (x,r) 


0 < C-V2ao)-^ < 


+ 


< 2(7^000, (x,r) e [xo,l] x [0,T*], 


r 2 (x,rl riix,rl 


(6.7) 

( 6 . 8 ) 

(6.9) 


r 2 (x,r) ri(x,r) 

(p.(p')-T'^ < < (p'lp.)-')'"'. e [0,x„] X [0,r.]. (6.10) 


with which, we can show the uniqueness of the solutions. Indeed, From fl3.4l) and using 
Young’s inequality, we have 


A. 

dr 


p^ ^R^dx = / (2p^ ^RRr + {6 — l)pi ^pirRjdx 
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= (1-6') / p?(rf ^ui)a:dx + 2 / p{ ^R{— - ^-^)dx 


r2 


U) 


<e / pl-^—dx + C, / pl-^Wdx, 


( 6 . 11 ) 


and 
d 

dr 


p{ ^Q^dx =2 / '’p(pir - P2T)dx + {9 -3) I Pi ^pirQ^dx 


X-3 


= (3-9) I Pi ^ni)a;p^(ia: + 2 / pi ^p(pi(rf - ^ 2 (^ 2 '' 


9-3 „2/N-l„ 




<e / Pi ^—dx + e / Pi'^Vf^ + / p® + Q / p® ^R^dx, 

Jo ri Jq Jq Jq 

( 6 . 12 ) 

where e > is chosen later and Cg > 0 a constant. 

From the equation fl3.4P o and boundary condition, we get 


^ f ^u'^dx= f {-6'pi+^(rf ^Mi)^(rf ^a;)^, - 6'p2+^(rf ^U2)x(r2 ^(^)x}dx 


+ / {pi'(rf-ia;).-pI(r2^-M4^^ 

Jo 

+ (A^-1) [ {pl(r^~‘^uicj)^ - pl(r^~‘^U2u)^}dx. 

Jo 

Using the similar argument as that in Lemma [3.21 and Cauchy-Schwartz inequality, we have 


d 1 
dr In 2 


u'^dx <e / p?’'’^ri^^ ^^u^dx + e / Pi ^^dx 


+ C, / pr P dx + Ce / Pi^R^dx, 

Jo Jo 


'0 


u 


(6,13) 


where £ > 0 small enough and Cg > 0 is a constant. 

Apply the Gronwall’s inequality to the summation of (I6.1ip - (l6.13p . we can hnally obtain 

f (vj^ + Pi~^R^ + Pi“^p^)(x, T)dx 

Jo 

<C [ (tc^ + pi“^i?^ + pi“^p^)(x, 0)dx = 0, (x,r) e [0,1] x [0,T*], 


which implies (pi,Mi,ri) = (p 2 ,U 2 ,r 2 ). 
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